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The polarization of prompt J/ip at the Tevatron is calculated within the nonrela- 
tivistic QCD (NRQCD) factorization framework. The contribution from radiative 
decays of P-wave charmonium states decreases, but does not eliminate, the trans- 
verse polarization at large transverse momentum. The prediction agrees with mea- 
surements from the CDF collaboration at intermediate values of pt, but disagrees 
at the large values of pt measured. 
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I The NRQCD factorization approach to inclusive quarkonium productionS 

lO ■ makes the remarkable prediction that in hadron collisionsL t / PC = quarko- 

nium states should be transversely polarized at large px □• In high-energy pp 
collisions, the dominant contribution to the charmonium production rate at 
large pr comes from gluon fragmentationa. At leading order in a s , a QQ 
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pair with small relative momentum created by the virtual gluon is in a color- 
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octet Si state with the same transverse polarization as the almost on-shell 
gluon. Due to the approximate heavy quark spin symmetry of NRQCD, 1 
quarkonium states from the pair should have a large transverse polarization 
at sufficiently large px- Recent measurements at the Tevatron by the CDF 
! collaboration seem to be in contradiction with this predictiona. 

A convenient measure of the polarization is a = (T— 2L)/(T+2L), where T 
and L are the transverse and longitudinal polarization fractions in the hadron 
■ CM frame. The variable a describes the angular distribution of leptons from 

the decay of the J /if) with respect to the J / ip momentum. Beneke and Roth- 
stein studied the dominantfragmentation mechanisms for producing longitudi- 
nally polarized 1 statestj. For charmonium production at the Tevatron, one 
should also take into account the fusion contributions from parton processes 
ij — » cc + k since fragmentation does not yet dominate for most of the pt 
region. The polarization variable a for direct J/ib and direct tp' at the Teva- 
tron have been predicted by Beneke and Kramer la and by LeibovichQ. They 
predicted that a should be small for pr below about 5 GeV, but then should 
begin to rise dramatically. On the contrary, the CDF data shows no sign of 
transverse polarization of direct rp 1 at large pr>3. 

The CDF collaboration has also measured the polarization of prompt J/ipa 
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(i.e. J/ip's that do not come from the decay of B hadrons). The advantage 
of this measurement is that the number of J /if) events is larger than for tp 1 
by a factor of about 100. On the other hand, theoretical predictions of the 
polarization of prompt J/ -0 are complicated since the prompt signal includes 
J/fkl$ that come from decays of the higher charmonium states Xcti Xc2, and 
if)' uC. The polarization of J/tp from ip' not via x.c is straightforward to cal- 
culate, since the spin is unchanged by the transition. The polarization of J/ip 
from Xc and of J/tp from ip' via Xc is more complicated, because the Xcj's 
are produced in various spin states and they decay into J/tp through radiative 
transitions ES. Therefore, it is interesting to investigate the cascade effect in 
this problem. 

The NRQCD factorization formula for the differential cross section for 
the inclusive production of a charmonium state H of momentum P and spin 
quantum number A has the schematic form 

da H ^ = da cS ^ (O^), (1) 

where the implied sum on n extends over all the color and angular momentum 
states of the cc pair. The cc cross sections da CCn , which are independent of 
H, can be calculated using perturbative QCD. All dependence on the state 
H is contained within the nonperturbative NRQCD matrix elements(ME's) 
(On X ^). The Lorentz indices, which are suppressed in (|l|), are contracted 
with those of da CCn to give a scalar cross section. The symmetries of NRQCD 
can be used to reduce the tensor ME's (On X ) to scalar ME's (O^) that are 
independent of P and A. One may calculate the cross section for polarized 
quarkonium once the relevant scalar ME's are known. 

In pp collisions, the parton processes that dominate the cc cross section 
depend on pt- If pt is of order m c , the fusion processes dominate. These in- 
clude the parton processes ij — > cc + k, with i, j = g,q,q and q = u,d, s. At pt 
much larger than m c , the parton cross sections are dominated by fragmentation 
processes. To get a consistent LO evaluation of polarization parameter a in 
this region, one should take care. Since we are interested in the LO prediction 
of a, we should first get a LO prediction of longitudinal fraction L/(T + L). 
The LO contributions to L are order-a^ processes g —> cc^r g, while the only 
contribution to T is the order-a s process g — > cc"8( 3 <Si) d. Inclusion of the 
longitudinal part in the fragmentation process is crucial in the calculation of 
a even though it has only a small contribution to the total cross section. The 
fragmentation function is evolved using the standard homogeneous timelike 
evolution equation. Since ccs( 3 Si) — » ipx(nS) is included in both the fragmen- 
tation and the fusion process, we interpolate between the fusion cross section 
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at low pr and the fragmentation cross section at high pT^'til- In all the other 
cc — > H channels, we use the fusion cross section. 

The color-singlet ME's can be determined .-uhenomenologically from the 
decay rates for tp(nS) — » £ + £~ and Xc2 — > 77^3. The color-octet ME's are 
phenomenologicaLfiarameters and they are determined from production data 
at the Tevatron EJtl. In case of Xc, we have a good agreement with the CDF 
data onihe-ratio Xd/Xc2&- More details on the analysis methods are explained 
in Refs.liJ'til. There are numerous theoretical uncertainties in the polarization 
calculation. In. our analysis, we allow the variations in PDF (CTEQ5L and 
MRST98LO) factorization and fragmentation scales (/it/2 — 2 fix), charm 
quark mass (1.45 — 1.55 GeV), respectively, where fix = (4m^+p|,) 1 / 2 . Wealso 
take into account the uncertainties in the ME's, (Os( 1 So)) and (Os( 3 -fo)) 

Let us compare our results for a with the CDF data. We present our 
result in the form of an error band obtained by combining in quadrature all the 
theoretical errors described previously. In Fig._l(a), we compare our result for 
direct ip' as a function of px with the CDF dataQ. But the error bars in the CDF 
data are too large to draw any definitive conclusions. We next consider the 
polarization variable a for prompt J/ip. The method fox.calculating cascade 
effect in prompt J/ip polarization is explained in Ref. 113. In Eig. 1(b), we 
compare our result for a as a function of px with the CDF datacl. Our result 
for a is small around px = 5 GeV, but it increases with px- Our result is in 
good agreement with the CDF measurement at intermediate values of px, but 
it disagrees in the highest px bin, where the CDF measurement is consistent 
with 0. The solid lines in Fig. 1(b) are the central curves of a for direct J/tp 
and for J/ip from Xc- The a for direct J/ip is smaller than that for direct tp' , 
because the ME's are differentlUJ. In the moderate px region, the contributions 
from and from Xc add to give an increase in the transverse polarization of 
prompt J/ip compared to direct J/ip. In the high px region, the contributions 
from ip' and Xc tend to cancel. 

The CDF measurement of the polarization of prompt J/ ip presents a seri- 
ous challenge to the NRQCD factorization formalism for inclusive quarkonium 
production. The qualitative prediction that a should increase at large px 
seems inescapable. However, it is still worthwhile to investigate the NLO ef- 
fect in the fragmentation processes, since the dramatic discrepancy appears at 
large px- Recently, Braaten and Lee performed the full NLO calculation of 
the colpr-octet 3 Si gluon fragmentation function for polarized heavy quarko- 
niumlix It is not sufficient to get NLO prediction of a, but one can first study 
the production rate at large px in NLO accuracy. 

In Run II of the Tevatron, the data sample for J/tp should be more than an 
order of magnitude larger than in Run I, which will allow both the production 
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Figure 1: Polarization variable a vs. for (a) direct j/i' and (b) prompt J/?/) compared to 
CDF data. 



rate and the polarization to be measured with higher precision and out to 
larger values of pr- If the result continues to disagree with the predictions 
of the NRQCD factorization approach, it would indicate a serious flaw in our 
understanding of inclusive charmonium production. The predictions of low- 
order perturbative QCD for the spin-dependence of cc cross sections could be 
wrong, or the use of NRQCD to understand the systematics of the formation 
of charmonium from the cc pair could be flawed, or m c could simply be too 
small to apply the factorization approach to the charmonium system. 
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